Journal of Cellular Biochemistry

Infiltration of the endothelial layer of the blood-brain barrier by leukocytes plays a critical role in health and disease. When passing through the endothelial layer during the diapedesis process lymphocytes can either follow a para-cellular route or a transcellular one. There is a debate whether these two processes constitute one mechanism, or they form two evolutionary distinct migration pathways. We used phylogenetic analysis, HH search, ancestor sequence reconstruction together with functional specificity and positive selection analysis to investigate this intriguing question further. We found that the two systems share several ancient components, such as RhoA protein that plays an important role in controlling actin movement in both mechanisms. However, some of the key components differ between these two transmigration processes. CAV1 genes emerged during Trichoplax adhaerens and it was only reported in trans-cellular process. Para-cellular process core proteins had at least two distinct starting points. First, during drosophila emergence, Tre1 which is homologous to melatonin GPCR receptor diverged. Secondly, PECAM1 emerged from FASL5/3 during elephant shark divergence. Lastly, both systems employ late divergent genes such as ICAM1 and PECAM1. Taken together our results suggest that these two systems constitute different yet interconnected mechanisms of immune cells infiltrations of the brain. Our analysis indicates that this system coevolved with immune cells, evolving to a higher level of complexity in association with the evolution of the adaptive immune system.

BCOR (BCL-6 corepressor)-rearranged small round cell sarcoma (BRS) is a rare soft tissue tumor, mostly featuring the BCOR::CCNB3 fusion, with other fusions like BCOR::MAML3, BCOR::CLGN, ZC3H7B::BCOR, KMT2D::BCOR, CIITA::BCOR, and RTL9-BCOR also reported. BCOR, a Polycomb Repressive Complex 1 (PRC1) component, influences histone modifications. It dimerizes with Polycomb group RING finger homolog (PCGF1) via its PCGF ubiquitin-like fold discriminator (PUFD) domain interacting with PCGF1s RING finger and WD40-associated ubiquitin-like (RAWUL) domain. We used various in silico tools to explore the impact of fusion events on BCORs functionality and RAWUL-PUFD dimer binding affinity. Changes were found in the domain landscapes, physicochemical properties, GO terms and significant increases in the disordered regions within the PUFD domain of the fusion proteins. Structural predictions indicated modified intermolecular contacts (ICs) and a significant reduction in binding affinity in fusion protein RAWUL-PUFD dimers. These findings align with expression data showing PRC1-regulated gene upregulation in BRS, likely due to reduced RAWUL-PUFD binding affinity, impacting dimer formation and PRC1 assembly. Our findings enhance the understanding of BRS oncogenesis and identify potential therapeutic targets.

This study presents a comprehensive bioinformatics analysis of the origin and structural complexity of thyroglobulin (TG). We examine the structural and evolutionary conservation of TG in Petromyzon marinus (sea lamprey) by reconstructing its complete TG sequence. Based on genomic data, we assembled a 2,831-amino-acid sequence (TGPM), identifying key TG domains and generating a homology-based PDB model. Additionally, we detected a second TG transcript in sea lamprey, designated TGPM1746 Comparative analysis across 38 representative vertebrate species--including mammals, birds, reptiles, amphibians, ray-finned fishes, and jawless vertebrates--reveals that all TG domains are conserved throughout vertebrate evolution. Despite substantial divergence in overall amino acid sequences, tyrosine residues and cysteines--both essential for TG function--remain highly conserved. TG emerges as a structurally complex, multidomain protein featuring a conserved disordered segment at its C-terminus. This region includes the terminal portion of the ChEL domain and the hormonogenic site responsible for triiodothyronine (T3) synthesis, likely contributing to the conformational flexibility required for hormone production. We further propose an evolutionary model in which a nidogen-like precursor--defined by the presence of TG type 1 modules--may have acted as the ancestral source of this essential repetitive motif within TG structure. Through genetic rearrangements and duplication events, a proto-TG likely arose, potentially shaped by environmental pressures such as ionizing radiation. Successive duplications expanded the TG architecture, culminating in the formation of 11 TG type 1 modules. The final evolutionary stage involved the integration of TG type 3 and TG type 2 modules, followed by the fusion of the ChEL domain, which enhanced TG secretion and thyroid hormone biosynthesis. Our findings demonstrate that the TG complexation process is fully established in lampreys and has remained remarkably conserved across vertebrate evolution.

Metaxin 3 genes are shown to be widely conserved in vertebrates, including mammals, birds, fish, amphibians, and reptiles. Metaxin 3 genes, however, are not found in invertebrates, plants, and bacteria. The predicted metaxin 3 proteins were identified by their homology to the metaxin 3 proteins encoded by zebrafish and Xenopus cDNAs. Further evidence that they are metaxin proteins was provided by the presence of GST_N_Metaxin, GST_C_Metaxin, and Tom37 protein domains, and the absence of other major domains. Alignment of human metaxin 3 and human metaxin 1 predicted amino acid sequences showed 45% identities, while human metaxin 2 had 23% identities. These results indicate that metaxin 3 is a distinct metaxin. A wide variety of vertebrate species--including human, zebrafish, Xenopus, dog, shark, elephant, panda, and platypus--had the same genes adjacent to the metaxin 3 gene. In particular, the thrombospondin 4 gene (THBS4) is next to the metaxin 3 gene (MTX3). By comparison, the thrombospondin 3 gene (THBS3) is next to the metaxin 1 gene (MTX1). Phylogenetic analysis showed that metaxin 3, metaxin 1, and metaxin 2 protein sequences formed separate clusters, but with all three metaxins being derived from a common ancestor. Alpha-helices dominate the predicted secondary structures of metaxin 3 proteins. Little beta-strand is present. The pattern of 9 helical segments is also found for metaxins 1 and 2.

The insulin superfamily proteins (ISPs), in particular, insulin, IGFs and relaxins are key modulators of animal physiology. They are known to have evolved from the same ancestral gene and have diverged into proteins with varied sequences and distinct functions, but maintain a similar structural architecture stabilized by highly conserved disulphide bridges. A recent surge of sequence data and the structures of these proteins prompted a need for a comprehensive analysis which connects the evolution of these sequences in the light of available functional and structural information and their interaction with cognate receptors. This study reveals a) unusually high sequence conservation of IGFs (>90%), which has never been reported before. In fact, it was interesting to observe that the functional domains (excluding signal peptide) of human, horse, pig and Ords kangaroo rat are 100% identical. (b) an updated definition of the signature motif of the relaxin family (c) a non-canonical C-peptide cleavage site in a few killifish insulin sequences and so on. We also provide a structure-based rationale for such conservation by introducing a concept called binding partners imposed evolutionary constraints. Furthermore, the high conservation of IGFs appears to represent a classic case of resistance to sequence diversity exerted by physiologically important interactions with multiple partners. Furthermore, we propose a probable mechanism for C-peptide cleavage in killifish insulin sequences.

The naked mole-rat (NMR, Heterocephalus glaber) is of significant interest to biogerontological research, rarely developing age-associated diseases, such as cancer. The transmembrane glycoprotein CD44 is upregulated in certain cancers and CD44 cleavage by a disintegrin and metalloproteinase 10 (ADAM10) regulates cellular migration. Here we provide evidence that mature ADAM10 is expressed in NMR primary skin fibroblasts (NPSF), and that ionomycin increases cell surface ADAM10 localization. However, we observed an absence of ADAM10 mediated CD44 cleavage, as well as shedding of exogenous and overexpressed betacellulin in NPSF, whereas in mouse primary skin fibroblasts (MPSF) ionomycin induced ADAM10-dependent cleavage of both CD44 and betacellulin. Overexpressing a hyperactive form of the Ca2+-dependent phospholipid scramblase ANO6 in NPSF increased phosphatidylserine (PS) externalization, which rescued the ADAM10 sheddase activity and promoted wound closure in NPSF in an ADAM10-dependent manner. These findings suggest that dysregulation of ADAM10 shedding activity is due to a deficient PS externalization in NMR.

SARS-CoV-2 is a novel virus causing mainly respiratory, but also gastrointestinal symptoms. Elucidating the molecular processes underlying SARS-CoV-2 infection, and how the genetic background of an individual is responsible for the variability in clinical presentation and severity of COVID-19 is essential in understanding this disease. Cell infection by the SARS-CoV-2 virus requires binding of its Spike (S) protein to the ACE2 cell surface protein and priming of the S by the serine protease TMPRSS2. One may expect that genetic variants leading to a defective TMPRSS2 protein can affect SARS-CoV-2 ability to infect cells. We used a range of bioinformatics methods to estimate the prevalence and pathogenicity of TMPRSS2 genetic variants in the human population, and assess whether TMPRSS2 and ACE2 are co-expressed in the intestine, similarly to what is observed in lungs. We generated a 3D structural model of the TMPRSS2 extracellular domain using the prediction server Phyre and studied 378 naturally-occurring TMPRSS2 variants reported in the GnomAD database. One common variant, p.V160M (rs12329760), is predicted damaging by both SIFT and PolyPhen2 and has a MAF of 0.25. Valine 160 is a highly conserved residue within the SRCS domain. The SRCS is found in proteins involved in host defence, such as CD5 and CD6, but its role in TMPRSS2 remains unknown. 84 rare variants (53 missense and 31 leading to a prematurely truncated protein, cumulative minor allele frequency (MAF) 7.34x10-4) cause structural destabilization and possibly protein misfolding, and are also predicted damaging by SIFT and PolyPhen2 prediction tools. Moreover, we extracted gene expression data from the human protein atlas and showed that both ACE2 and TMPRSS2 are expressed in the small intestine, duodenum and colon, as well as the kidneys and gallbladder. The implications of our study are that: i. TMPRSS2 variants, in particular p.V160M with a MAF of 0.25, should be investigated as a possible marker of disease severity and prognosis in COVID-19 and ii. in vitro validation of the co-expression of TMPRSS2 and ACE2 in gastro-intestinal is warranted.

Breast Cancer Metastasis Suppressor 1 (BRMS1) expression is associated with longer patient survival in multiple cancer types. Understanding BRMS1 functionality will provide insights into both mechanism of action and will enhance potential therapeutic development. In this study, we confirmed that the C-terminus of BRMS1 is critical for metastasis suppression and hypothesized that critical protein interactions in this region would explain its function. Phosphorylation status at S237 regulates BRMS1 protein interactions related to a variety of biological processes, phenotypes [cell cycle (e.g., CDKN2A), DNA repair (e.g., BRCA1)], and metastasis [(e.g., TCF2 and POLE2)]. Presence of S237 also directly decreased MDA-MB-231 breast carcinoma migration in vitro and metastases in vivo. The results add significantly to our understanding of how BRMS1 interactions with Sin3/HDAC complexes regulate metastasis and expand insights into BRMS1s molecular role, as they demonstrate that BRMS1 C-terminus involvement in distinct direct protein-protein interactions.

M. tuberculosis is one of the most successful human pathogens causing tuberculosis that leads to highest daily morbidity worldwide. The evasion of the host immune responses is an important strategy that M. tuberculosis adopts. MprA (Rv0981), the response regulator of two component system is known for DNA binding activity in the pathogen and its role in persistent infection in the host. MprA is recognized as a late stage antigen during infection. A variant form of the protein MprA with G70S polymorphism (MprA*) is observed in one of our local and in several global clinical isolates of M. tuberculosis. Here we report the nuclear localization of MprA and MprA* in differentiated macrophages. MprA and MprA* increase the expression of TGF-{beta} and IL-10, the immune suppressive cytokines in THP-1 derived macrophage cells. Concurrently the phago-lysosome fusion is significantly reduced as shown by infection with M.bovis BCG. We show that single nucleotide variation in clinical isolates lead to quantitative variations resulting in host immune suppression and support the survival and persistence of the pathogen.

The Spike glycoprotein receptor-binding domain (RBD) of SARS-CoV-2 mediates the viral particles binding to the angiotensin-converting enzyme 2 (ACE2) receptor on the surface of human cells. Therefore, Spike-ACE2 interaction is a crucial determining factor for viral infectivity. A new phylogenetic group of SARS-CoV-2 (lineage B.1.1.7) has been recently identified in the COVID-19 Genomics UK Consortium dataset, which features an amino acid substitution in the Spike RBD (N501Y mutation). Infections with the SARS-CoV-2 lineage B.1.1.7 have been overgrowing in recent weeks in the United Kingdom, indicating an even greater spread capacity than that seen with previous strains of the novel coronavirus. We hypothesized that this rapid spreading/infectivity of the B.1.1.7 lineage might be due to changes in the interaction force between the mutant Spike RBD and ACE2. This study employed in silico methods involving mutagenesis (N501Y mutation) and interface analysis focusing on the Spike RDB-ACE2 interaction. The results showed that the SARS-CoV-2 N501Y mutant (lineage B.1.1.7) establishes a more significant number of interactions relating to the mutant residue Y501 (Spike RDB) with residues Y41 and K353 (ACE2). This finding shows that the increased infectivity of SARS-CoV-2 lineage B.1.1.7 is associated with the interaction force between the Spike RBD Y501 mutant residue with the ACE2 receptor, which in this strain is increased.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) utilizes TMPRSS2 receptor to enter target human cells and subsequently causes coronavirus disease 19 (COVID-19). TMPRSS2 belongs to the type II serine proteases of subfamily TMPRSS, which is characterized by the presence of the serine-protease domain. TMPRSS4 is another TMPRSS member, which has a domain architecture similar to TMPRSS2. TMPRSS2 and TMPRSS4 have been shown to be involved in SARS-CoV-2 infection. However, their normal physiological roles have not been explored in detail. In this study, we analyzed the amino acid sequences and predicted 3D structures of TMPRSS2 and TMPRSS4 to understand their functional aspects at the protein domain level. Our results suggest that these proteins are likely to have common functions based on their conserved domain organization. Furthermore, we show that the predicted 3D structure of their serine protease domain has significant similarity to that of plasminogen which dissolves blood clot, and of other blood coagulation related proteins. Additionally, molecular docking analyses of inhibitors of four blood coagulation and anticoagulation factors show the same high specificity to TMPRSS2 and TMPRSS4 3D structures. Hence, our observations are consistent with the blood coagulopathy observed in COVID-19 patients and their predicted functions based on the sequence and structural analyses offer avenues to understand better and explore therapeutic approaches for this disease.

Tunneling nanotubes (TNTs) are thin cytoplasmic protrusions involved in long-distance cellular communication. The presence of TNTs has been found in vivo and in vitro studies in non-small cell lung cancer (NSCLC). Cancer cells transport a range of organelles and signalling molecules along TNTs, to confer a survival phenotype for the recipient cell, contributing toward chemoresistance and malignancy. Despite its important role in cancer progression, the molecular mechanisms underlying TNT formation is not well defined. Within the tumour microenvironment (TME) of NSCLC, hepatocyte growth factor (HGF) and its receptor, c-Met, are mutationally upregulated causing growth, and invasion. In this study, we report a novel crosstalk between HGF/c-Met and {beta}1-integrin involved in the formation of functional TNTs in A549 cells. Through pharmacological inhibitor studies, we discovered Arp2/3 complex, MAPK and PI3K pathways were activated downstream of this crosstalk signalling axis. Furthermore, paxillin was recruited during this key process, localising at the protrusion site of HGF-induced TNTs, and therefore serving as the central link between the upstream and downstream regulators involved. Overall, these results demonstrate a novel strategy to inhibit TNT formation in NSCLC through targeting the HGF/c-Met and {beta}1-integrin signalling axis, thus highlighting the importance of personalised multi-drug targeting in NSCLC.

Host-pathogen interaction is the best example of an evolutionary arms race where pathogen and host continuously coevolve to survive and exert negative effects on each other. The adaptability of both host and pathogen is critical for this association. In this study, we explored the association of severe acute respiratory syndrome (SARS) coronaviruses (CoVs) with their human host from the genomic and evolutionary perspectives based on a comparative analysis of SARS and MERS coronaviruses. We observed that human proteins that are part of the SARS-CoV2-human association are enriched in hubs and bottlenecks. Again, these proteins take part in more protein complexes and show slower evolutionary rates compared to the human proteins associated with the two other coronaviruses, SARS-CoV and MERS-CoV. Moreover, the human proteins involved in the interaction with SARS-CoV2 are mostly longer proteins harboring long intrinsically disordered stretches and a higher level of disordered protein binding sites. Codon usage analysis revealed that the novel coronavirus is least adapted to codons used in housekeeping and lung-specific genes, compared to the other two coronaviruses. We conclude that the signatures showed by the SARS-CoV2-human protein interaction network revealed the viruss association with vital human proteins and pathways, via interactions mediated by protein complexes and intrinsically disordered protein binding sites, which may have assisted the higher infectivity of SARS-CoV2 in its human host than the other two less-virulent human coronaviruses, despite having a lower optimization to its hosts codons.

The P4-ATPase family of phospholipid flippases plays a critical role in the maintenance of membrane asymmetry and consequently, various roles in cellular protein traffic and eukaryotic homeostasis. Currently, several structures of these (usually heterodimeric) phospholipid flippases have been resolved, along with extensive biochemical characterization of the substrate transport properties. However, an essential subfamily of monomeric phospholipid flippases, the P4B-ATPases, remains to be characterized in depth. While these P4B-ATPases appear to have similar lipid transport properties to their heterodimeric counterparts, the P4A-ATPases, the basis of their existence as monomers is currently unknown. In addition, the unique features of this group have yet to be comprehensively analyzed since the resolution of one P4B-ATPase member. In this study, we investigated the divergence of P4B-ATPases from other P-type ATPases using a structure-based sequence analysis approach. This led to identification of features unique to P4B-ATPases, as well as previously unidentified conserved properties of P4-ATPases. The results of this study provide a basis for further studies on P4-ATPases to characterize conserved properties of this group that supersedes substrate specificities. MANUSCRIPT INFORMATION- Number of manuscript pages: 20 (double-spaced). - Number of Figures: 5. - Number of Tables: 3. - Supplementary materials include 2 supplementary tables.

Non-coding, single-stranded RNA molecules known as microRNAs (miRNAs) regulate gene expression via mRNA degradation after transcription. As a result, they affect a number of pathways in organisms that are important for both health and disease. miRNAs can be utilized as potential diagnostic, prognostic, and therapeutic biomarkers for neurodegenerative diseases such as Amyotrophic Lateral Sclerosis (ALS). Neuronal cells are highly dependent on mitochondria, and mitochondrial dysfunction has been linked to neurodegenerative diseases. Pathological changes in ALS are associated with disruptions in mitochondrial structure, bioenergetics, and calcium homeostasis. In this study, we used an in silico approach to identify miRNAs associated with mitochondrial dysfunction in ALS based on target genes that are implied in both ALS and mitochondrial dysfunction. A literature search revealed the genes SOD1, FUS, TARDBP, C9orf72, CHCHD10, OPTN, VCP, TBK1 and BCL2 that cause mitochondrial dysfunction and are involved in the pathogenesis of ALS. Pathway enrichment analyses using Enrichr, g:Profiler, and CROssBAR tools confirmed that the identified genes have significant associations with ALS, mitochondrial dysfunction, and neuron differentiation. In silico miRNA predictions have been made using the databases miRWalk, miRTargetLink, TargetScan, and miRNet. A Venn diagram tool was used to select common miRNAs, and finally 28 miRNAs were discovered. One set of 28 miRNAs were subjected to set analysis using the miRNet and TAM tools for functional and enrichment analyses, respectively. In both databases, three common miRNAs, hsa-miR-9-5p, hsa-miR-141-3p and hsa-miR-125b, were found to be linked to ALS.

The monoamine oxidase (MAO) gene family encodes for enzymes that perform the oxidative deamination of monoamines, a process required to degrade norepinephrine, serotonin, dopamine, and other amines. While mammalian MAO enzymes, MAO A and MAO B, have been extensively studied, the molecular properties of the other family members are only partly uncovered. This study aims to explore the evolution of monoamine oxidases, emphasizing understanding the MAO gene repertoire among vertebrates. Our analyses show that the duplication that gave rise to MAO A and MAO B occurred in the ancestor of tetrapods, between 408 and 352 million years ago. Non-tetrapod jawed vertebrates possess the ancestral preduplicative condition of MAO A/B. Our results also identified a new family member, MAO C, in non-tetrapod jawed vertebrates. Thus, most jawed vertebrates possess a repertoire of two MAO genes, MAO A and MAO B in tetrapods and MAO A/B and MAO C in non-tetrapod jawed vertebrates, representing different MAO gene lineages. Jawless vertebrates possess the ancestral condition of a single copy gene, MAO A/B/C. Enzymatic assays conducted on the MAO recombinant enzymes of the Indo-Pacific tarpon show that both proteins, MAO A/B and MAO C, have enzymatic and molecular properties more similar to human MAO A, with the former featuring a strikingly higher activity rate when compared to all other MAO enzymes. Our analyses underscore the importance of scanning the tree of life for new gene lineages to understand phenotypic diversity and gain detailed insights into their function.

SARS-CoV-2 is highly contagious and can cause acute respiratory distress syndrome (ARDS) and multiple organ failure that are largely attributed to the cytokine storm. The surface coronavirus spike (S) glycoprotein is considered as a key factor in host specificity because it mediates infection by receptor-recognition and membrane fusion. Here, the analysis of SARS-CoV-2 S protein revealed two B56-binding LxxIxE-like motifs in S1 and S2 subunits that could recruit the host protein phosphatase 2A (PP2A). The motif in S1 subunit is absent in SARS-CoV and MERS-CoV. Phosphatases and kinases are major players in the regulation of pro-inflammatory responses during pathogenic infections. Moreover, studies have shown that viruses target PP2A in order to manipulate hosts antiviral responses. Recent researches have indicated that SARS-CoV-2 is involved in sustained host inflammation. Therefore, by controlling acute inflammation, it is possible to eliminate its dangerous effects on the host. Among efforts to fight COVID-19, the interaction between LxxIxE-like motif and the PP2A-B56-binding pocket could be a target for the discovery and/or development of a bioactive ligand inhibitor for therapeutic purposes. Indeed, a small molecule called Artepillin C (ArtC), a main compound in Brazilian honeybee green propolis, mimics the side chains of LxxLxE motif. Importantly, ArtC is known, among other effects, to have anti-inflammatory activity that makes it an excellent candidate for future clinical trials in COVID-19 patients.

The lipolytic enzymes of Mycobacterium tuberculosis play a critical role in immunomodulation and virulence. Among these proteins, PE11 which also belongs to the PE/PPE family, is the smallest ([~]10.8 kDa) and play a significant role in cell wall remodelling and virulence. PE11 is established to be an esterase, but its enzymatic and structural properties are not yet characterized. In this study, using homology modelling we deduced the putative structure which shows the presence of both -helix and {beta}-sheet structures which is in close agreement with that observed by CD spectra of the purified protein. PE11 was found to contain a GX3SX4G motif homologous to canonical GxSxG motif present in many serin hydrolases. The catalytic triad appears to be located within this motif as substitution of Serine26 and Glycine31 residues abrogated its enzymatic activity. Gel-filtration chromatography data indicate that PE11 possibly exists as dimer and tetramer showing positive cooperativity for binding its substrates. In addition, PE11 esterase activity was found to be critical for cell wall remodelling, antibiotic resistance and conferring survival advantages to M. tuberculosis. Our data suggest that PE11 can be targeted for designing potential therapeutic strategies.

Severe respiratory syndrome coronavirus 2 (SARS-CoV-2) infection presents an immense global health problem. Spike (S) protein of coronavirus is the primary determinant of its entry into the host as it consists of both receptor binding and fusion domain. While tissue tropism, host range, and pathogenesis of coronavirus are primarily controlled by the interaction of S protein with the cell receptor, it is possible that proteolytic activation of S protein by host cell proteases also plays a decisive role. The host-cell proteases have shown to be involved in the proteolysis of S protein and cleaving it into two functional subunits, S1 and S2, during the maturation process. In the present study, the interaction of S protein of SARS-CoV-2 with different host proteases like furin, cathepsin B, and plasmin has been analyzed. Incorporation of the furin cleavage site (R-R-A-R) in the S protein in SARS-CoV-2 has been studied by mutating the individual amino acid. Our results suggest the polytropic nature of the S protein of SARS-CoV-2. Our analysis indicated that a single amino acid substitution in the polybasic cleavage site of S protein perturb the binding of cellular proteases. This mutation study might help to generate an attenuated SARS-CoV-2. Besides, targeting of host proteases by inhibitors may result in a practical approach to stop the cellular spread of SARS-CoV-2 and to develop its antiviral.

The SARS-CoV-2 outbreak was recently declared a worldwide pandemic. Infection triggers the respiratory tract disease COVID-19, which is accompanied by serious changes of clinical biomarkers such as hemoglobin and interleukins. The same parameters are altered during hemolysis, which is characterized by an increase in labile heme. We present two approaches that aim at analyzing a potential link between available heme and COVID-19 pathogenesis. Four COVID-19 related proteins, i.e. the host cell proteins ACE2 and TMPRSS2 as well as the viral protein 7a and S protein, were identified as potential heme binders. We also performed a detailed analysis of the common pathways induced by heme and SARS-CoV-2 by superimposition of knowledge graphs covering heme biology and COVID-19 pathophysiology. Herein, focus was laid on inflammatory pathways, and distinct biomarkers as the linking elements. Finally, the results substantially improve our understanding of COVID-19 infections and disease progression of patients with different clinical backgrounds and expand the diagnostic and treatment options.